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THE  DISTRIBUTION  OF  A  HEAVY  POLYDIbPERSED 
AEROSOL  IN  A  TURBULENT  ATMOSPHERE  AT  A  LONG  DISTANCE 
FROM  AN  INSTANTANEOUS  POINT  SOURCE 

[Following;  is  a  translation  of  an  article  by 
I.  L»  Karol'  and  A.  Ya,  Pressman  in  Inaheaerno- 
FizichesHiy  Zhurnal  (Engineering-Physical  Jour¬ 
nal),  Vol.  II,  No,  9,  Minsk,:  1959,  pages  83“ 

91.  J 

Thia  work  analyses  certain  qualitative  laws  of  the  distribution 
in  th®  atmosphere, and  of  th©  gravitational  deposition  to  the  sur¬ 
face  of  the  earth, of  an  aerosol  that  Is  heterogeneoua  with  respect 
to  its  particle  composition  and  that  la  at  a  considerable  distance 
from  its  source.  It  dlsoueaesi  the  case  of  an  instantaneous  point 
source  that  is  located  at  heights  h  above  the  surface  of  the  earth} 
in  additlosri,  h  Is  mush  greater  than  the  thickness  of  the  surface 
layer  of  the  atmosphere. 

1,  Within  the  limits  of  the  semi-empirical  theory  of 

turbulent  diffusion,  the  non-stationary  process  of  dispersion 

of  a  jHo  nodiaper  si  ve  aerosol  in  the  field  of  a  Hnwlni? 

Is,  as  is  known  fl],  described  by  the  solution  of  tne  iollov^lng 

equation  with  constants 


_L  - fjc  »^S!i 

di  .  dx  dz 


«? 


-r  K 


-f 


W 


which  satisfy  the  initial 

=»Q5(4S(y)5{« -ft) 

and  boundary  conditions 

npH  X*  4*  y®  ~i"  -►00, 


(1) 

(2) 

(3) 


1 


where  z  --  o  Is  t-he  surface  of  the  earth?  qw  is  the  yQl.um3  eoa- 
centration  of  a  monodispar sive  aerosol 

w.  rL  the  constant  velocity  of  a  wind  f 

x\r>  nx*  Kv  K~t  and  K?  are  constant  coefficients  of  turDalent 

diffusion  of^  iha  aerosol  in  horizontal  and 

tions*”  S(x)  is  the  delta^funct ion  of  Dirac;  Q  is  the  total 

LinSSt’of  aeroLl  given  off  by  the  source;  /S  is  the  factor  oi 

proportionality,  between  the  flow 

7-0  and  the  vo3.ume  cone  entrap  ion  ol  the  feerusol,  «nx^n 

Jharaoterlze  the  Interaction  of  particles  of  the  aerosol  with 

the  3“J„f,^„^S*1,iitlc.n  (1)  Which  satlsflas  the  con- 

ditiona  of  (2)»  (3)  and  (4)  has  tne  form  llj. 


Sir^tf^V/CJCj^K 


X 


X  « 


w.f 


(s+fi)‘ 


2e 


h  4^2  -\~  «>  t 

“2vt;F, 


(5) 


where 


m 


1/^1^  f  w 

2p -  ®  > 0  ;  ^  (C)  =  ~54— e-^erfcC 


C  r"’'"  d  t, 


As  in  reference  [2],  let  us  assume  that  the  jerosol, 

heterogenous  with  respect  to  composition,  maj  0. 

described  by  incomplete  ?/rlT  “ihoc^e  deF 

respect  to  the  velocities  of  precipitation  w  L3J,  wnose  deu 

sity  of  distribution  has  the  forms 

«&>  ”1  rt 

(6) 


Niw)  = 


W 


'■m 


m 


where  parameter  Wm  characterizes 

fraction  prevailing  with  rart3» 

n  is  the  degree  of  homogeneity  of  the  distribation  of 

cles  of  the  aeroBol  with  respect  wO  the  Yelocitl  s  ,  *  (  } 

#K — “"If — anTliT'niaT'Tse’Tonsld©^  constant,  if  the  source  is 
loSatS  moh  SighL  than  the  thickness  of  the  surface  layer 

si  Thj  SSaitloris.  0  la  sufficient  for  the  unlquenase  of  the 
soltition  of  the  problems  (1)  -  (4), 


the  gamma  function. 

In  distribution  (6)  let  us  note  two  extreme  oases; 

a]  the  aerosol  is  homogenous,  i.e.,  it  consists  of  parti^ 
which  nreolnitate  in  the  atmosphere  with  the  same  tsio 
cliy  w  ~  wj  in  this  case,  n  infinity  and  (6)  transforms 

Into 

A'  (?»)  S  iw  —  Wg,) ; 


b)  virtually  non-precipitating  particles  prevail  in  the 
zeroBol/ ~  Q-,  in  this  case,  assumlrig  Wjn  =  nwg  and 
transf orxttlng  in  t^)  tov/ards  the  limit  Wi. wh  n  0  we  oOtci>xn. 


1  ^ 

N{'W  — ~e  «'■'« 


(8) 


In  eonneotion  with  the  heterogeneity  of  the  aerosol, 
it  is  necessary  to  take  into  account  the  dependence  Oi  the 
value  of  |5  and- the  coefficients  of  turbulent  diffusion  on 
the  velocities  of  fall  of  particles  w. 

It  is  natural  to  assume  that  the  condition  of  absorp-^ 
tion  is  realized  for  comparatively  heavy  particles  on  uiie  sun 
face  of  the  earth,  i.e.,  generally  speaking 


on  the  other  hand,  at  a  great  dlstance^  from  the  source 
^ with  respect  to  the  wind),  when  the  practically  weightless 
coBiPonent  of  the  aerosol  Is  dispei’sed,  the  interaction  oi  uh© 
flow  of  the  aerosol  v/ith  the  surface  of  the  earth  (  z  =  0) 
will  be  determined  by  the  value; 

\w^0 


¥• 


shall  examine  two  cases; 


i )  —  C)»  i«e», 


.  the  "welf^btless"  component  of  the  aero- 

sol  is  repelled  by  the  surface  of  the  earth;. 

?) /ie  >  0,  l.e.,  on  the  surface  of  the  earth  partial  or 
total*' =  <so  )  absorption  of  the  "welghtlexSB  conponerit  of 
the  aerosol  takes  place. 

The  coefficients  of  turbulent  diffusion  Kx,  y,  z 
the  sa.’.ne  way  as  in  [2],  we  will  consider  to  ba  aepen- 
dent  on  w  in  such  a  waj!'  that 


K.t,  y.  z  (w)  -^  0  K  /G.  z  («)->•  K^.  y,  *  (0). 


-J 


Let  us  give  the  coefficients  Kx,  y,  z  smaller  values,  insofar 
as  we  regard  the  dispersion  of  ihe  virtually  non»prec iplta- 
tinpj  conmonent  of  the  aerosol,  wltn  large  values  of  time  t» 

2®  Let  us  hegin  to  obtain  formulas  for  the  volurce  and 
surface 'concentration  of  a  polydispersed  aerosol  at  a  large 
distance  from  the  source;  furthermore,  let  us  begin  from  the 
case  =  Ou  .Assumlna:  here  the  arbitrary  power  dependence 
/5(w)  :r  bw®,  from  conditions  ^  {0}  -“0  and^^J  »  0  we  will  find* 
that 

p  =--buK  (9) 

with  b  =  i  the  boundary  condition  (4)  has  the  appearance; 

==0  (10) 


dz 


b-O 


and  ai'‘xnifi8s  the  absence  of  a  turb'ale.nt  flow  of  the  aerosox 
through  the  surface  z  0.  This  boundary  condition  serves  in 
problems  (1)  -  (4)  for  considering  the  influence  of  the  sur¬ 
face  layer  of  the' atmosphere  on  the  distribution  of  the  set¬ 
tled  admixture  on  the  surface  of  the  earth,  ^ 

AssuminiS?  in  (5)  ~  ~  »  lot  u.s  xirrite  tne.i.ol,i.ox^iiig 

expression  'for '"the  volume  concentration  of  a  polydispersed  ae. 
rosol; 

q  {x,  )f,  z,  t)  ~  I  q^,  N  (w)  dw  = 

0 


2n. 


w 


V 


F. 

t 


« •‘V'  n 


Qexp  ^  — 


4KJ 


4K,t 


8 («  +  !)/*/» I 


X 


«2 


where 


X  1  I ,,  ■  ^  X  XF  -f  X)  { i  X. 

0 


(H) 


>  --  ®  /  t  .  , 

2}/  7c';  ' 


h 


2VKJ 

n  .  rw: 


h~-  ^  ^ 

Yv^i' 


K'‘m  i 


t 


(12) 


5h4ctuarry,""let  /S 


bwaTT  foA  fullfillment  under  all  condi' 


tioiis 


a;  =  2  bw®  -  w>0  it  is  necessary  that  m  -  1 


4. 


into 
l\,ie  s 

pandy 


where 


is  thi 

(11) 


whei^e 


first 


Let  us  construct  the  expansion  of  the  function  (ll)-- 
L  series  with  respect  to  powers  of  small  (for  large  va*« 

)f  t)  paraiaetersL,f/  and  ;  preliminarily,  let  us  ex~ 
with  respect  to  powers  or  the  function  f  % 

iiS) 


'(^(c®-f-x)= 

. .,  J.™ 

'  i  (k) 

5  function  of  Hermlte  CS];  after  which  let  us  transcrihe 
In  the  form 


\  V  t  j  I 


Qe 


X 


£f  -"’tl 


e  'ie 


^fnAh) 


4^ 


X 


*-i 


(I4» 


mf  >,•?  m 

4.o(y  -  f  =  2";™  f 

5  S  « 


s 

X«-»-»//,{X)/f_*(X)r-*^i^X. 


O 


:f'  «i« 


(15) 


iimm  * 

Limiting  ourselves  in  'both  obtained  series  to  only  the 
terms,  we  haver 

/«.« (e,)  ^  I X*  d  X  +  0(S,)  ^  -J  r  j  +  0  (It), 


/«+f,0  (^i)"  j  '  //_t  (X)e-'*£;x  +  0(U 


(16) 


.5. 


L1+.1  j 


4/1-4- 


ir 

2  \  2  / 


+  0(?i). 


Then  the  first  term  of  tha  ey^sansion,  of  the  function  (3.4) 
with  respect  to  powers  of  v/xXl.  assume  the  loi'm  of 


q(jc,  y,  X 


exp 


X 


4l'(x  ^  4/^ 


/-  ,  (ft  —  zf  n{h  —  z),  n*  j 

}  ^  ^  ^  iAji-wn  ^  ■  lunn.iin.iJrt  iftrviY"  1 

V'  '■■  wj  <rl.  y 


g  «  r  +’i 


X 


n 


-wr^e  Oil,) _ 


(17) 


Let  us  note  that  term-by-term  Integration  of  the  ae¬ 
rie's  for  4^  In.  (3  3)  is  justified,  since,  the  series  obtaxnea, 
by  substitution  of  (13)  into  (14)  agrees  uniformly  with  res¬ 
pect  to  A  in  X  3s  0.  ~  XU- 

In  the  case  under  considerat ion^  the  formuia  lor 
surface  concentration  of  an  aerosol  wliich  has  settled  on  (..ne 
earth  is  obtained  from  expression  (5)  by  the  follov/ing  means: 


™  r  ^  _5^_ -j- I  Af(?Kf)rfa'«=  i  ^{w)q^ 

J  I  J  i-o  x! 

0  t’’ 

/  2n  irjr-Y*' _  ^  . _ 


*«»o 


-{n+SV? 


X 


X 


/m-..  0  iw-  2  51  {“  t/v"iW"'44-i  #0 

j&mi 


dt, 


where 


4/4^  4/0  mj 


(IS) 


(19) 


and  Iv  firCi)  Is  determined  by  expression  (15).  Since  we  are 
analyzi-y  the  dispersion  of  the  aerosol  at  a  great  disbance 
‘^'Vrsm  ti-ir-}  sourco*  v/e  ar©  intCirBst^d  Iti  thB  yeIub  o  ix,  j)  xor 
iaJSe  values  of! L.  Therefore,  having  utiliz.ed  ^ae  toowr^Lj- 
Dlace  method  [5  I.  lot  us  construct  the  asympuotio  expanuiyO. 
of  the  integral  (18)  for  L  Infinity.  In  this  case,  as  lo 


6 


known  C53»  the  c oof i^lc lent s  of  the  asyaptotic  expansion  (18) 
are  determined  from  the  expansion  of  the  extraexponential 
factor  of  the  Integrand  expression  in  (18)  with^respect  to 
the  negative  powers  of  te  Thus.,  considering:  il6),  we  can 
write 

®  (Xy  ™  I  . . .  V  I  Q  “  i  i\  ""  *X. 

\  ‘^m  / 


2  k 


X 


“  _L  jbV 

ft.,/ 


.3,r(l±i) 


VK.H 


4*5 


t  *  (1 + di 


\  / 


Q 


Sit  («+3)  j  mrl) 


rt+S 

"T" 


X 


X  |ft{n+a>/'2(«  VWKr)  +  O  f  ^  4)/2 


(20) 


of  Macdonald  [4]®  Utilising  the 


where  K^jCz)  is  the  function, 
asymptotic  expansion  ILt(z)  for  large  values  of  argument  z,  we 
find' the  final  exp;resslon  for  the  sain  terra  of  the  asymptotic 


y) 


i-^vrT' 


" 

exp  - 

XU 

2K,-  “■* 

tVLjKr 

Y7Ky{n+Z)r. 

l  +  of  ^ 

1^) 

- 

« 

X 


i21) 


Let  us  exami.ne  the  case  of  partial  absorption  of 
the  "weightless"  component  of  the  aerosol  by  the  surface  of 
the  earth,  i.e*,  the  case  /S  (0)  >-0,  Here^  for  large  val^aes 
of  time  tj,  the  argument  of  the  function  i|/  in  expression  (b) 
Increases  limitleasly;  therefore,  for  construction  of  the  a- 
symptotlc  expansion  q(x,y,  z,t)  and  i7""(x,  y)  at  large  distan¬ 
ces  from  the  source,  let  us  utilise  the  asymptotic  expansion 
of  the  funetioB,  4'' (  fo:r  ^  ->  infinity  [4-3,  and  let  us  pr0“ 
sent  the  extrae:xponential  part  of  the  third  addend  of  expres¬ 
sion  (5)  in  the- form; 


«)# 


(22) 


where 


«,=2: 


(2^  -  1  YMl  4-  ife)!  {2K,Y  (fe  +  g)*  . 
2ife!  (/ —  fe)!  «)*+' 


«6 


A  -J-  2  .  2/Cj, 

I ;  flj  *=  — - 1 - ;  ' 

u>  w* 


ComeutiriK  the  volume  concentration  of  a  Poly^ispersed 
substitute  this  asymptotic  expression  for  t-^ 
initv  into  thA  third  addend  of  the  integrand  expression 
Ul)  and  a?ter  term-by-term  Integration  with  respect  to  w  we 

shall  obtain* 


%/  W  W  WA* 

/-F-y-t  '  exp  {-(X  -  utri4K>i  -  f  l±l^ 

qix,%zJ)’^\^—y  -^)  Q  8(r/)V.VX^ 


X 


(23) 


iasrl 


where 


d, 


r(rt4' 


namely 


+  l)x 

x2  (-2J^)’. 


(24) 


rTerm-bY^reTi^intgsr  in“"r22 )  with  respect  to  w  from  0  to 

tnfiSlt(  iflustified,  since,  (25)  by  virtue  of  f 

>0,  this  expansion  is  uniform  with  recpect  to  w  in  w  ^  0. 
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where  the  coefficients  b 
Sion 


a/  exp 


ar^e  determined  from  the  expan-"" 

<25) 


'h—z _ 

.  2Kt 

arastO 

w-O 


Limiting  ourselveSj  as  in  the  previous  case,  by  the 
main  terms  in  the  obtained  expansions,  we  vrill  find  from  ex¬ 
pression  (2i)  the  main  term  of  the  asymptotic  expansion  q(x, 
y,z,t)  With  t  -3»  infinity:* 


q{x,y,zJ) 


n  Q/C,"/=>exp  {-  jx  -  tttfjAKJ  -  yV^Kj]  ^ 


X 


(■ 


4  V  _  e 


m 


The  asymptotic  expression  for  the  surface  concentration 
(r(x,  y)  at  large  distances  from  the  source  Is  found,  in  the 
same  way  in  paragraph  2,  and  its  main  term  has  the  form  of 


a(x.y)-f~-y _ 

I  wj  8«V/<^yr(/j/2-f  I) 

m 

I 


(-t) 


X 


X  ^ 


t 


2 


['+"(w)] 


dt. 


Carrying  out  the  integration,  and  utilizing  the  asymptotic  re¬ 
presentation  of  the  Macdonald  function  Kjji(z),  we  have,  analo¬ 
gously  to  (21), 


/  n  \«-n  Q^"^'^exp[xif/2/Cr  — 

\  j  4  /k  r  {rt/2  +  1 )  (4/C 


(27j 


*In  this  case,  the  equality 


W,(  01  =  /ir2’' 


1 


is  utilized. 


9. 


Let  us  note  that  if  condition  of  absorption  of 

the  aerosol  by  the  earth  s  surface  is,  satisfied. 

0, 


nJ  = 

«  •  ♦ 


i  e  in  the  specific  casie  of  condition  (4)  ith/?  - /6d( 

aSSm^  'Sa  (?3)  is  eaisrtlllly  h^ilhed;  and'tae 

appearance  of. 


(x.y) 


QHKrg 


■n/a 


n  X'’^' 

— )  X 

®m/ 


<28> 


*«»o 


C  «=  (u^^fe  -  2raK*)/'K?*,  V'K,- 

4.  Wt  ue  analyse  the  for  lar- 

Froiu  formulas  <27),  it 

ge  values  of  the  value  ^  surface  concentra- 

source  with  respect  to  win  ,  k  decreases  with  an  Increase 
tlon  of  the  fallen  aerosol  for  x infi- 

nf  TT  S  r*  Or*<i  iOR  to  tho  6Xp02l^Btft2l.l  4.1  “hp  P-* 

llty  to  (21 W  (87),  the  factors  which  contain  x  will  he  a- 

symptotically  equal  to*. 


where 


4Kxl 


1  + 


4Kx 

(  , 

n%\ 

X* 

1  4/f, 

/ 

«  x*fi 

exp|x«/2/r^~-Kl/I7Cl 

0X'p  I  /2  Kx  ('-/ 


(29) 


«tbe  main  term  iVHr5Sy^p^S,g^,|;^?f 
passing  to  the  limit  tor Infinity. 
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Thus,  in  case  of  repulsion  of  the  "welshtless"  compur 
nent  of  an  aerosol  by  the  surface  of  the  earth  =  0.  fomu- 
la  (21)),  far  from  the  source  cT  (x,  y)  decreases  as 
and,  in  case  of  even  partial  absorption  >  0,  fomula  (27))^ 
decreases  as  Let  us  note  that  for  n  ^  infinity,  — 

=  Wai  0  and  A>  0,  i.e,,  in  case  of  precipitation  of  a  heavy 
homogenous  aerosol  on  the  absorbing  (even  if  partially)  sur¬ 
face  z  ~  0,  or"  (x,  y)  decreases  with  an  increase  of  x  according 
to  the  exponential  law  (as  expression: 


A'"*  exp 


A-"*exp  — 


for  small  w™),  and  not  according  to  the  square  law,  as  is  er¬ 
roneously  stated  In  reference  (6).  The  surface  concentration 
of  a  monodlspersed  aerosol  decreases  with  an  increase  oi  x  as 
x-2  only  In  the  case  of  Wjjj  ~  0  and  (partial  absorption 

of  the  "weightless”  aerosol)*  ,  .  ,  .  * 

As  is  seen  from  formulas  (21),  (27)  and  expressions 
(29),  for  X  infinity  the  main  terms  of  the  asymptotic 
panslons  (x,  y)  do  not  depend  on  Kx; 

fusion  In  the  direction  of  the  wind  does  not  inliuence  the 
distribution  of  the  surface  concentration  of  an  aerosol  pre¬ 
cipitated  at  large  distance  from  the  source,  under  any  boun¬ 
dary  conditions  on  the  surface  of  the  earth. 

The  main  terms  of  the  asymptotic  expansions  I x,  y;  in 

(21)  and  (27)  depend  essentially  on  the  value  of  K^;  i.e,,  the 
role  of  tur-bulent  diffusion  with  respect  to  the  vertical  is 
essential  in  preclplta.tlon  of  an  aerosol  at  large  distances 
from  the  source  in  the  direction  of  the  wind,  except  for  the 
case/3<,>0,  n  =  0  (distribution  of  the  type  (8)),  when  with 
wm  nwo  and  n  -»  0  the  main  term  in  the  right  side  (2f)  beco¬ 
mes  Independent  of  Kj-*  .  ,  x  . 

Thus,  some  quantitative  and  qualitative  characteristics 

of  dispersion  of  a  polycllspersed  aerosol  in  a  turbulent  medium 
may  differ  considerably  from  the  appropriate  regularltxes  ol 
the  distribution  of  a  homogenous  aerosol,  even  if  it  is  assumea 
(as  is  done  in  the  present  work,  as  well  as  In  12J>  tnat  va¬ 
rious  fractions  of  the  polydlspersed  aerosol  disperse  indepen¬ 
dently  of  each  other. 
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SUMMARY 

The  paper  deals  with  qualitative  and  quantitative  characteristics  o 
the  process  of  dispersion  of  polydisperse  aerosol  a  long  distance  downwind 
of  the  source. 

The  relationship  between  the  expressions  for  volume  concentration 
and  surface  concentration  of  polydisperse  aerosol  deposited  from  the 
atmosphere  ind  boundary  conditions  on  the  surface  of  the  earth  Is  deri¬ 
ved  in  pp.  1  3.  Thus,  formulae  for  volume  and  surface  concentrations 

are  shown  to  be  essentially  different  depending  on  whether  or  not  the 
.weightless*  component  of  polydisperse  aerosol  Is  reflected  by  the  sur¬ 
face  of  the  earth  (Formulae  17,  20,  21  and  26,  27,  28). 

A  detailed  qualitative  analysis  of  the  formulae  obtained  is  given  inf 
B-i  Ik  llli' the  or^'igisaX  .doouaient]* 
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